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Abstract. We compute the connected four-point correlation function of the 
primordial curvature perturbation generated during inflation with standard kinetic 
terms, where the correlation is established via exchange of a graviton between two 
pairs of scalar fluctuations. Any such correlation yields a contribution to the scalar 
trispectrum of the order of the tensor to scalar ratio r. This contribution is numerically 
one order of magnitude larger than the one previously calculated on the basis of scalar 
perturbations interacting at a point and satisfies a simple relation in the limit where 
the momentum of the graviton which is exchanged becomes much smaller than the 
external momenta. We conclude that the total non-linearity parameter generated by 
single- field models of slow-roll inflation is at maximum [tnlI ~ f- 
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1. Introduction 

Since the release of the first-year Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite data, the availability of high-resolution all-sky maps of the Cosmic Microwave 
Background (CMB) anisotropy has revolutionized cosmology as a quantitative science. 
Together with large galaxy surveys, the ever-improving accuracy of the microwave 
background data has yielded a great deal of precise information concerning conditions 
in the very early universe. Over the next several years a new generation of CMB 
instruments will return yet more data, which will fix the parameters of the standard 
model with even greater precision and may point the way to the discovery of new physics. 

Despite its role as the dominant carrier of information from the earliest times, the 
statistical properties of the CMB temperature anisotropy are remarkably simple: it is 
known to be statistically Gaussian to very high precision, greater than 0.1%. In recent 
years, however, the possibility that small non-Gaussian features might be present in 
these fluctuations has received increasing interest from both cosmologists and particle 
physicists. This is because primordial non-Gaussianities are one of our best probes 
of whatever complex physics was operating in the very early universe. In particular, 
non-Gaussian features discriminate sharply between the minimal implementation of 
inflation — driven by a single scalar field rolling down its potential under the influence 
of Hubble friction, known as "slow-roll" dynamics — and more complicated variants 
involving extra fields or dynamics of greater subtlety. 

Which observables are relevant in the study of non-Gaussianity? An exactly 
Gaussian process is one in which all statistical information is encoded in the two-point 
expectation value. Any observation of new information in the three- or higher ra-point 
correlation functions therefore constitutes an observation of non-Gaussianity. Other 
statistical measures exist, such as Minkowski functionals [Ij. However, at the present 
day the best confrontation of theory with observation is by means of the three-point 
correlation function of the curvature perturbation, (. 

To make this quantitative, one defines the bispectrum 5^ by 

(C(ki)C(k2)C(k3)) = (27r)35(^kj5^(ki,k2,k3) , (1.1) 

a 

where is a function of the triangle formed by the three momenta {ki,k2,k3}, and 
we have introduced a notation which will be used throughout this paper, in which 
Latin indices such as a,b,c, . . . range from 1 to the number of momenta carried by the 
correlation function in question. These are to be distinguished from Latin indices such 
as i,j, ranging from 1 to 3 which label the spatial dimensions. It is customary to 
introduce a dimensionless quantity, /nl, which generally depends on the momenta ka, 
and parametrizes the strength of the non-Gaussian signal. We define 

Bc^lhLY.Pdka)Pdh), (1.2) 

a<b 

where P^ki) is the power spectrum of C, such that (C(ki)^(k2)) = (27r)'^5(^^ ka)P^(fci). 
Occasionally, /nl is referred to as the non-linearity parameter ^ i3j. An /nl that is 
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independent of the momenta corresponds to a local type of non-Gaussianity, which arises 
from the real space relation 

c = c. + If^i^X - iO) ' (1-3) 

where (g is Gaussian. In this case, the bispectrum is maximized in the limit of one 
of the three momenta going to zero P]. This form of non-Gaussianity is produced by 
nonlinear gravitational evolution subsequent to horizon crossing [H [5] . 

On the basis of non-linear perturbation theory we might naively expect a non- 
linearity parameter of order unity, i.e., /nl ~ 0(1). However (neglecting post-processing 
in the later universe from the ubiquitous non-linearities of gravity [6l [71 [H [9l [10]) in 
single-field slow-roll inflation /nl is very small — in fact, it is of the same order as the 
slow-roll parameters [HI [12], which are required to be much smaller than unity. Such 
a small non-linearity is beyond our present experimental reach, which is |/nl| ^ 100. 
However, in the near future the Planck satellite may be sensitive to /nl of order unity 
[13] . Any detection of a primordial | /nl | larger than order unity will strongly disfavour 
minimal inflation. 

It is possible that interesting information might be encoded in the higher-order n- 
point functions, such as the connected part of the four-point function [H]. Constraining 
such an expectation value from CMB data is extremely challenging and at the time of 
writing no such analysis has been completed for the WMAP data, although a study on 
small angular scales has been carried out for the BOOMERanG data [15] and an all-sky 
analysis of the four-point function exists for the older COBE satellite [16] . However, as 
computing power increases and algorithms become more efficient, it is quite plausible 
that it will be possible to extract four-point correlations for the WMAP data, and that 
this will be possible for Planck. In view of future constraints from CMB experiments 
and elsewhere (see for instance Refs. [IZl [IE]), the purpose of this paper is to clarify 
the predictions of single-field slow-roll models of infiation for the four-point correlation 
function of the curvature perturbation (. 

The four-point correlation function is conventionally parametrized in terms of the 
trispectrum, written T^, and defined by [T9l [20] 

(C(ki)C(k2)C(k3)C(k4)) = (2vr)35( J] k,)T^(ki, k^, kg, k4) , (1.4) 

a 

where by convention only the connected part of the expectation value on the left-hand 
side is considered. The disconnected part is present even in an exactly Gaussian model, 
and forms a background from which the genuinely non-linear signal must be extracted. 
Unlike the bispectrum, the trispectrum typically depends on the relative orientations of 
the ka in addition to their magnitudes. If observation forces us to abandon the minimal 
scenario of infiation, then the trispectrum is a ready-made observable which could be 
consulted for guidance concerning the modifications which are required. 
The trispectrum can be parametrized as [21] 

54 

= TNL Yl Pdkab)Pdh)Pdkc) + —9NL Yl Pd^a) Pcih) Pd^c) , (1-5) 

b<c a<b<c 
a^b,c 
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Figure 1. Possible planar momentum quadrilaterals. Note, however, that there is no 
need for the quadrilateral to lie in a plane, and in general it is a three-dimensional 
object. In (a), there is no particular relationship among the sides of the quadrilateral. 
This is the general case. In (b), the side associated with k4 is taken to zero length, 
causing the quadrilateral to degenerate to a triangle. In (c), adjacent sides possess 
equal magnitudes and their directions are becoming opposite, a limit we refer to as the 
folded kite. In (d), opposite sides possess equal magnitudes, yielding a parallelogram. 
The shapes described in (c) and (d), are therefore precisely dual to each other, and 
describe the same physics. 



where kat = \i^ab\ and kab = + kf,. This definition is motivated by the fact that, 
for tnl and (7nl which are independent of the momenta, this is the most general local 
parametrization, corresponding to a non-Gaussianity of the form 

c = c, + ^ {r^i^'f" (c| - (c,^)) + |^7j.r'a (1-6) 

in real space. As for the bispectrum, this simple momentum dependence characterized 
by constant r^l^^ and g^^l^^ is generically produced by gravitational evolution subsequent 
to horizon crossing. 

More generally, the classical gravitational evolution generates a curvature 
perturbation which can be expanded in powers of (g with definite momentum- 
independent coefficients. In such a case, if the perturbations are generated by only 
one field, one expects that [21] 

^local / ^ f local \ /i ^\ 

^NL = 1 5/NL I ■ 

However, the curvature perturbation during infiation is subject to quantum interference 
effects around the time of horizon crossing, and does not admit a simple parameterization 
in terms of momentum independent-coefficients [22l H] . It follows that in general there 
need be no relationship between tnl and /nl- 

There are two interesting degenerate configurations where a non-Gaussianity of 
local form maximizes the trispectrum. The first is the so-called "squeezed" limit, where 
one of the four momenta is taken to zero. In this limit, the quadrilateral formed by the 
momentum vectors k^ degenerates into a triangle. (See Fig. [Tl^a)-(b).) For instance, by 
taking ^4 -C ^1,^2, k^, one obtains 

Tl""^^' = (^2r^T' + ^9'^l'^ Pdh)m^^ . (1.8) 

Thus, T^{f^^ and (7n£^' contribute equally to the trispectrum in this limit. 
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The second interesting limit occurs when the magnitude of the sum of two momenta 
is taken to zero, so that kij = |kj + kj| ^0 for some i ^ j- We refer to this hmit as 
the counter- collinear case because for each momentum there is another one with equal 
magnitude and opposite direction. One can visualize the quadrilateral formed by the 
momentum vectors in this limit in two different but equivalent ways. If we choose to 
order the momenta so that the counter-collinear pairs are adjacent, then we generate 
what can be called a "folded kite," shown in Fig. [TJ^c). On the other hand, if we 
choose to order the momenta so that the counter-collinear pairs are on opposite sides 
of the quadrilateral, then it becomes planar and yields a parallelogram, as in Fig. [D^d). 
Changing the order of the momenta does not affect the trispectrum, which is a function 
only of the intrinsic geometry of the k^; for the trispectrum, this means it depends on 
a total of three magnitudes and three angles. Therefore, the planar shapes given in 
Fig. [Hc)-(d) are precisely dual to each other. The local trispectrum simplifies in the 
counter-collinear limit. By taking ki2 -C fci ~ A;2, ^3 ~ A;4 one obtains 

T^°^-i(ki, k2, k3, k4) = AT'^l^'P^{h2)Pdki)Pdh) ■ (1.9) 

Thus, only t^l'^^ contributes to the trispectrum in this limit. 

As discussed above, it does not always happen that a model of the early universe 
predicts a curvature perturbation whose non-linearities are precisely of the local form 
described by Eqs. (11.31) and (11.61) . To accommodate this situation, one can take tnl 
and s^NL as momentum dependent; however, in this case the separation between tnl 
and (7nl becomes ambiguous. Indeed, several authors [23], [5], [2l] have used a momentum 
dependent tnl alone to parametrize the full trispectrum. The division of the trispectrum 
into terms proportional to tnl and (^nl is helpful where these parameters can be made 
approximately momentum- independent, but this is not always the case when the non- 
Gaussianity is dominated by quantum interference effects around the time of horizon 
exit, as in single-field inflation. For this reason we will mostly avoid the use of tnl 
or ^fNL in this paper, giving expressions for the trispectrum directly. However, as we 
will see, the contribution of graviton exchange to the trispectrum is dominant in the 
parallelogram limit, in such a way that we can describe it as a momentum independent 

_local 
^NL • 

In perturbation theory our expectation would again be tnl ~ fi'NL = ^l}-)- However, 
it has been shown that, in slow-roll inflation with canonical kinetic termj§| tnl exhibits 
a similar suppression by one power of the slow-roll parameter e = ~H/H^, giving 
Tnl ^ r/50 [28], where r is the ratio of tensor to scalar amplitudes, r = 8Pt/P(, 

§ The calculation of the trispectrum in single-field models with non-canonical kinetic terms was first 
performed, neglecting metric perturbations, in Ref. [26j and later, more generally, in Ref. [27j . In these 
models tnl can be of order unity or larger. The contribution from the graviton exchange is suppressed 
by slow-roll with respect to the one due to the contact interaction and it has been neglected in those 
references. Note that the calculation performed in this paper, i.e., of the trispectrum from slow-roll 
inflation models based on a field with standard kinetic term, requires one to pursue the calculation to 
one order higher in slow-roll than was done in Refs. [26l[2I]- Thus, those references do not include our 
final results. 
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Figure 2. In (a), exchange of a graviton (represented by a wavy line) leads to 
correlations among four scalar fluctuations (represented by straight lines). In the limit 
where the graviton which is exchanged becomes extremely soft, so that |ki + ^ 0, 
one can think of such correlations as being mediated by fluctuations on top of a modified 
background which carries a classical gravitational wave. The same interpretation holds 
in (b), where the exchange is mediated by a scalar particle. 



and in single-field inflation r = 16e. It follows that the non-linearity parameter tnl is 
roughly of the same size as /nl- (In practice, of course, the trispectrum is much harder 
to detect than the bispectrum because an extra power of the power spectrum is 
equivalent to suppression by ~ 10^^°.) However, one may have some reservations about 
this conclusion. In particular, the calculation reported in Ref. |28] ignored nonlinear 
interactions mediated by scalar or tensor fluctuations. In terms of diagrams, this is 
equivalent to neglecting the two diagrams where the correlation among the four scalar 
modes is due to the exchange of scalar or tensor quanta. The latter diagram, which 
turns out to be more important, is shown in Fig. ([2]) (a), whereas the scalar exchange is 
described in Fig. El^b). 

The vertices of these nonlinear interaction can be computed straightforwardly using 
the action for Einstein gravity coupled to a single scalar fleld, which was studied to 
third order by Maldacena [TT]. For the vertex which describes a contact interaction 
between three quanta of the curvature perturbation, (, the appropriate term in the 
Lagrangian is typically suppressed by at least one power of the slow-roll parameters, 
i.e., ~ 0(e, 77) ■ and therefore 



where 77 is the second slow-roll parameter, deflned as 77 = MpV" /V, where we have 
introduced the reduced Planck mass Mp = (SvrG)"^/^. Thus, bearing in mind that a 
four-point correlation established by scalar exchange must contain two vertices of this 
sort, we expect that the contribution to the trispectrum from this process is of the order 
of two powers of slow-roll, i.e., tnl ~ 0{e'^,ri'^,eri). This is suppressed by one power 
of slow-roll compared with the scalar contact interaction and can therefore be ignored. 
However, things are different when we consider the exchange of a graviton. In this case 
the interaction Lagrangian between a graviton and two scalars is not suppressed by any 
powers of slow- roll parameters at all, i.e., C^ccl ^CC ~ -^7^^ therefore 



The power spectrum of gravitational waves is smaller than that of ^ by a factor of r, and 
therefore this contribution is of order r. It follows that the effect of graviton exchange is 






1/2 
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comparable to the contact contribution studied in Ref. |2S]. This observation was first 
made by Arroja & Koyama [STjljjl 

In the present paper, we will calculate the contribution to the four-point function 
of scalar curvature perturbations which arises from correlations mediated by graviton 
exchange. We explicitly carry out this calculation for single-field models, but our 
result is easily generalized to many fields and our final results will apply to multi- 
field examples. Like many previous calculations of the non-linearities among scalar 
fluctuations at horizon crossing, this calculation is somewhat involved. Fortunately, 
in the counter-collinear limit {i.e., parallelogram or folded-kite), the amplitude and 
momentum dependence of the contribution to the trispectrum from graviton exchange 
can be computed straightforwardly by a simple semiclassical argument, which provides 
an important check of our computation. 

In the next section we explicitly compute the contribution to the trispectrum 
from graviton exchange, using the uniform curvature gauge, while in ^ we discuss 
the counter-collinear limit and check our calculation. In §1] we give an expression for 
the total trispectrum, including the contribution of Ref. [28] and in ^we conclude with 
a discussion. Finally, in Appendix A we report the calculation of a long double integral 
used in ^ 

Throughout this paper, we use fundamental units in which Planck's constant H and 
the speed of light c are set to unity, so h = c = 1. We choose to set Mp ^ = SnG equal 
to unity. Our metric convention is (— , +, +, +). 



2. Trispectrum from graviton exchange 

In this section we want to evaluate the contribution to the four-point function of 
curvature perturbations from the exchange of a graviton. This is the process shown 
in Fig. El^a), which involves a third-order interaction among scalar fluctuations and 
tensor perturbations of the metric. The description of this interaction requires that we 
expand the action to third order in such fluctuations. For simplicity, we will restrict to 
a single scalar field only. The calculation can easily be extended to more than one field, 
and indeed our final result for the trispectrum of the curvature perturbation ( remains 
valid in the multi-field case with flat field space metric. 

We begin with Einstein gravity coupled to a scalar field, which gives a total action 

S=\j V^[R- {V<P)' - 2V{<j>)] , (2.1) 

II In Ref. 27] it was also claimed that the discussion of Ref. [29] neglected tensor modes. However, the 
tensor contribution originally derived by Maldacena [11], which was discussed again in Ref. [27], was 
in fact included in Ref. 29J, where it was shown to confirm their conclusions. 

One might also worry that the graviton-scalar-scalar vertex could dominate corrections to the scalar 
power spectrum, which arise via loops which exchange virtual quanta. The correction to the scalar 
power spectrum due to loops from the graviton-scalar-scalar interaction has been computed in [30j . 
where it was shown that it is of the same magnitude as that coming from scalar self-interactions 

[SD [Sana EH [35]. 
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where R is the spacetime Ricci scalar and V{(j)) is an arbitrary potentiaL Using the 
Arnowitt-Deser-Misner method (or so-called "ADM formalism") [36], one can derive 
the action for coupled scalar and perturbations from Eq. (12. ip to any given order using 
an iterative procedure. The ADM line element is given by 

ds"^ = -N'^dt^ + a^{t)hij{dx' + N'dt){dx^ + N^dt) , (2.2) 

where a{t) is the scale factor, and A^, A^* are the lapse and the shift functions. The three- 
dimensional metric hij encodes scalar and tensor fluctuations in the spatial geometry. 
When the ADM metric ansatz is inserted into Eq. (12. ip . the lapse and the shift functions 
grange multipliers: the field equations obtained by minimizing the action in A^ 
and A^* give the constraint equations which follow from the Einstein equation, while (f) 
and hij are the dynamical degrees of freedom. 

As in Ref. |28], we choose to work in the uniform curvature gauge. In this gauge 
we take the scalar field and the spatial metric to be parametrized by 

= 0(t) + (^(t,x) , hij = a{t){e^)ij . (2.3) 

Here V5(t, x) is the inflaton field fluctuation, and is a traceless and divergenceless 
tensor fluctuation, obeying di'~^ij = •ju = 0. Our eventual aim is to compute the 
four-point function of the comoving curvature perturbation, TZ, which on super-Hubble 
scales is equal to the curvature perturbation on uniform density slices, ( [371 EH 139] . 
However, the calculation is considerably simplified if we choose to calculate the four- 
point function of if at horizon crossing and then relate it to the four-point function of 
the uniform density curvature perturbation (. Indeed, as shown in Ref. [H], in virtue of 
the smallness of the slow-roll parameters at Hubble crossing, it is possible to reproduce 
the correct late-time behaviour of ( using the free-field solution in de Sitter spacetime 
for ip. At late times ( is constant and, at lowest order in slow-roll, it can simply be 
normalized to ip at horizon crossing. 

The free action, which describes the evolution of independent scalar and tensor 
fluctuations, is 



So = d^x dt 



■2 



djjpdip\ 1 f . . _dkjijdkjij\ 



:2.4) 



where we are using the convention that a repeated spatial index in the lowered position, 
such as ^ or j, denotes summation with the flat spatial metric bij. The free fields derived 
from this action satisfy the equations of motion obtained by varying Lp and 7jj, and can 
therefore be built out of mode functions Uk(t), U^it) and 7|(t), 7|*(t) which solve the 
classical equations of motion, 

Ukiv) = ^il-tkr])e'''\ (2.5) 



/2A;3 

= ^{1 - tkv)e'''' (anys), (2.6) 



where rj is the conformal time, defined by dr] = dt/a, and is the Hubble parameter 
evaluated at Hubble crossing for the mode k; an asterisk "*" as superscript denotes 
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complex conjugation. The index s labels the two possible polarization states "+" 
and "x" of a spin two gravitational excitation. Introducing creation and annihilation 
operators a^, for the scalar quanta we obtain 

^[ak^.(t) + al,f/,*(t)]e^'^-, (2.7) 

whereas, with corresponding operators b^'^, b^'^ for each polarization of the gravitational 
fluctuation, we find 

[K'^UWkit) + &+k4(-k)7r(i)] e^'-^ • (2.8) 

The polarization tensors efj are chosen to satisfy the transversality and tracelessness 
conditions e|j(k) = k^e^j{k) = 0, together with a completeness relation obtained by 
tracing over spatial indices, e|j(k)e*j^'(k) = 25ss' 

We have now reduced the problem to calculating the four-point function of 
fluctuations cp in the inflaton field at horizon crossing, where the correlations are 
produced via the graviton exchange diagram. To proceed, we need only the relevant 
third order term in the action, which describes the graviton-scalar-scalar vertex in 
uniform curvature gauge. This is pT] . 

Sz = \j d^x dri a'^-f'^d^ipdjip . (2.9) 

The correct approach to computing expectation values in an expanding or curved 
background spacetime is the so-called in-in or Schwinger-Keldysh formalism p0 ll^H2] . 
In this formalism the expectation value of any operator O (evaluated at time r]^) is given 
by 

(0|t(c»(?/,) e-*/-*oc'^'?[-f^i('^'^+'T+)--f^i(<^''/'"'T")]| |0) 

{n\ oiv.) \n) = ^ r . ' (2-io) 

(0| T < e~* ''f^^'^''^ '"^ \ |0) 



where \Q) is the vacuum of the interacting theory while |0) is the vacuum of the free 
theory, T is the time ordering operator and Hi is the interaction Hamiltonian given, in 
our case, b}0 

Hi = -^j £x a'^Y^diipdjip . (2.11) 

The matrix element (12.101) describes a system evolved from its initial vacuum state at 
rj — s> — oo to f] = rj^ with an operator inserted at rj, and back again from 77* to —00, 
with a set of "+" fields on the increasing-time contour and a set of "— " fields on the 

% In general, the interaction Hamiltonian does not have to be given by the sign reverse of the interaction 
Lagrangian. This is only true in a limited class of theories without derivative interactions, which does 
not include the case of coupled matter and metric fluctuations during inflation. It follows that extra 
terms can be generated in the transformation to the interaction Hamiltonian [431 125j . For models with 
non-minimal kinetic terms, this can lead to non-negligible contributions to the trispectrum, as discussed 
in v3 of Ref. [26]. In the case of slow-roll inflation with canonical kinetic terms, however, these effects 
can be ignored. 
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decreasing-time contour. Eventually, 77* will be taken to and the final result is finite 
in this limit. 

The contractions between different pairs of the two types of field perturbations yield 
four kinds of propagator, 

x') ^ (0| T [^^{x)^^ix')] |0) . (2.12) 

It is convenient to rewrite this expression in Fourier space, where the propagator takes 
the form 

G^Hx,x') = I (0e^"(— ')G^±(r^,V) , (2.13) 
and can be written as 

G,-iv, i) = GUv, v'Mv' -v) + GUv, v'Mv - V') , 

G,^r^,v')=G>{v,v'), 

Gt-iv,v')=GUv,v') ■ (2.14) 
The ">" and "<" functions are defined by 

G>{v,v')=U,{vmiv'), 

G<{v,v')=U;{v)U,{v') . (2.15) 
Similarly, for the tensor modes we have a propagator F^^^^^ 

Fgr^Hx,x') ^ (0| Tbti^hti^')] |0) , (2.16) 
which goes over to the Fourier space representation 



:pss'±±( I 
ijlm y-^j-^ 



__e^k.(x-x'),.^.(k)e^^(_k)F,^±±(r^,V) . (2.17) 



As before, explicit expressions can be given in terms of ">" and "<" functions, leading 
to 

Fr-(r/,r/') = (r^, V)0(V - v) + V)0(^ - v') , 

Fr+(r^,V) = i^r(^,^'), 

Fr-(r/,r/') = Fr(r],V) , (2.18) 

where one has 

Friri,v')=ir{vH{v') ■ (2.19) 

Expanding the exponential in Eq. fl2.10p to second order, using the interaction 
Hamiltonian given by Eq. (12. lip , and making all possible Wick contractions among the 
fields gives the contribution to the four-point function from graviton exchange (which 
we label with the superscript "GE"). Since there are two vertices in the diagram, the 
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result depends on two time integrals. Using Eqs. fl2.12p and fl2.16p to carry out the Wick 
contractions, together with the Fourier space definitions (12.131) and (12.171) . leads to 



12 



X G++(r^„ r^)G++(r/., v')Gl^{ri.. v') 

- Gl+iv*, v)Gt;{v*, v)Gt{v., v')Gl-iv*,v')F^^,-iv, v') 

- G+-(r/., v)Gt{v*, v)Gt;iv*, r/')^47(^' 

+ 23 permutations . (2.20) 

Using G^^{7]^,7]) = G^{r]^,r]) for r]^: > r], F^^{7],r]') = F^{r],r]') for r] > r]', G> = G^* 
and other equivalent relations, this can be simplified to give 

(y^k.y^k.y^kaV'kjr = (2vr)353(^kJ^e|^.(ki2)A;ifc^2<™(k34)4fcr T § /' ^ 

as ~°° 

X Im [G> (r/., r/)G,> (^*, ^)] " Im [G.^l^*, (r?*, r/')i^.t ^0] 

+ 23 permutations . (2.21) 

To proceed, one makes use of Eqs. (12.151) and (12.191) for the propagators G^{rj^,ri') 
and F^^{r], rj') in this equation, together with the free mode functions for the scalar field 
and tensor fiuctuations, Eqs. (12. 5p and (12. 6p . In total, for the four-point function we 
obtain0 

m 



(V'k.y^k.^ka^k,)?^ = (27r)^^(^k,) ; ■ 5^e|,.(ki2)6L(k34)A:iF24^r -^1234 



+ 23 permutations 
where X1234 is defined as 

2^1234 = 



(2.22) 



■00 J ^00 v'"^ 



X Im 



1 + iA;3?7*)(l + ik47]^){l - ik^7]'){l - ik^r]') 



(2.23) 



X (1 + ikuv){'^ - zA;i2V)e'^'^'"'^*^^'''+''*^e-*(''~^')^'i 

and kt is the sum of the magnitudes of the momenta, kt = ^a- 

Let us work out the permutations in more detail. Since e^j is a symmetric tensor 
and Eq. (12.220 is invariant under the exchanges 1 2 and 3 4, one can rewrite this 
equation as 

AH^ 1 

(¥'ki¥^k2</'k3¥'k4)™ = (27r)=^(5(^k,) * ■ W^U^-^)^ln^i^cd)''aKk^K 



-^abcd 



s a<b 
c<d 



(2.24) 



+ An equivalent — and computationally equivalent — interaction picture formalism was derived by 
Weinberg [43]. With this formalism the trispectrum can be computed straightforwardly using Eq. (1) 
of this Ref., expanded to second order in the interactions. 
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where the sum is only over different indices a, b, c, d, and we have collected an overall 
symmetry factor of 4 which takes into account the exchanges a ^ b and c ^ d. 
Rewriting this summation somewhat more explicitly, one arrives at the formula 

1 



E 



pr4(''l2)C(k34)*i*^*3*r 



4 



(^1234 + 2^3412) 



(^1324 + ^2413) 



(2^1423 + 



2314J 



(2.25) 



The calculation of X1234 is extremely lengthy but straightforward, and is given 
explicitly in Appendix A There are divergent terms in the limit r]^ together with 
a logarithmic dependence on the momenta, but these peculiarities cancel between X1234 
and X3412 to give an answer which is finite at late times and depends rationally on the 



X1234 + X3412 — 

kik2 



ki + ^2 



^-34 



(^34 + ki2){al^ - 2634) + kl^{k:i + k^) 



+ 



&34 , , ki2 

ki2^ 

034 ai2 



k3k4 



ku^]{^ + - 



0-34 



ai2 



+ (1,2^3,4) 



+ (1,2^3,4) 



k 



12 



0-120-34:kt 

where we have used ki2 



612634 + 2kU n k,) f 1 + + jr^] 1 , (2.26) 

V V«^t ^12^34 ktai2a3ij _ 



kt 012^^34 
/C34 and we have defined 



Oiab — ka -\- k[) -\- k, 



ah 



Jab 



(2.27) 



{ka ~\~ k^jkab ~\~ kakf) . 

To simplify Eq. ( 12. 25^ further, let us work out the polarization factors. To start, 
we will express '^s^ijO^^'i)^tm0^3'i)kikik^3k2' in terms of the relative angles between the 
ka and ki2. The polarization tensors can be rewritten as 

(2.28) 



p.p. 

where e and e are orthogonal unit vectors perpendicular to ki2 . It is convenient to write 
the ka in a spherical coordinate system having {e, e, ki2} as basis (we have denoted 
k = k//c). In this coordinate system one obtains 

ka = fca(sin 9a cos 0a, sin da sin 0a, cos 9a) , (2.29) 



where cos = ka ■ ki2 and cos ( 



k„ ■ e. Thus 



eljk\kl = kik2 sin 9i sin 6^2 cos(0i + ^2) 



(2.30) 

€ijk\k2 = kik2 sin 9i sin 6*2 sin(0i + ^2) , (2-31) 

with an analogous relation holding for e^jk^kl and e^jk^kl. Since the projections of 
ki and k2 (k3 and k4) on the plane orthogonal to ki2 have the same amplitude but 
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opposite directions, then sin 6^2 = /i:isin6'i and 02 = 0i + tt (fc4sin^4 = k^smO^ and 
04 = 03 + vr). Using these relations, one finally finds 

^ e|^.(ki2)eL(k34)A;ifc^4fcr = klkl sin^ 9, sin^ ^3 cos 2xi2,34 , (2.32) 

s 

where Xi2,34 = 0i — 03 is the angle between the projections of ki and k3 on the plane 
orthogonal to ki2. This is also the angle between the two planes formed by {ki, k2} and 
{k3, k4}. This expression can easily be generalized to the polarization factors in the last 
two lines of Eq. ([225]). 

We now wish to convert our prediction for the scalar four-point correlation function 
due to graviton exchange into a prediction for the correlation function of the curvature 
perturbation, (. In order to do so, it is sufficient to use the linear relation between ( 
and the scalar field perturbation ip at Hubble crossing, i.e., ( = (p^,/V2e. As shown 
in Refs. [251 HU [21], higher order terms in the relation between ( and generate 
contributions to the trispectrum that are suppressed by slow-roll parameters with respect 
to the contribution from graviton exchange. Thus, using Eq. f l2.32|) and analogous 
expressions for the other polarization factors, we can rewrite Eq. (12.251) to give 



k^k^ - - - - 

-^[1 - (ki ■ ki2)'][l - (k3 ■ ki2)'] cos 2X12,34 " (2:1234 +23412) 
^^12 

+ T3^[l - • kl3)'][l - (k2 ■ ki3)2] COS 2X13,24 " (^1324 + X2413) 
^^13 



2u2 



+ ^[1 - (ki ■ ki4)^][l - (k2 ■ ki4)^] COS 2X14,23 " (Xl423 + X2314) 
fcl4 



(2.33) 



3. The counter-collinear limit 



The calculation of the four-point correlation induced by the graviton exchange process 
is rather complicated, and the final expression given in Eq. (12.331) — combined with 
Eq. (I2.26|) — is not simple to interpret. One would therefore like an independent check 
that the result is correct. It is possible to achieve this in the counter-collinear limit, 
where k^b = k^ + k^ — and we can attach a simple physical meaning to the exchange 
process: in this limit, the four-point function simply expresses the correlation between a 
pair of two-point functions, for instance (CkiCk2) and (Ck3Ck4) which is induced by a low 
frequency gravitational wave of momentum ku <^ /ci ~ k2, ks ~ k^. This gravitational 
wave crossed the Hubble radius much earlier than any of the ka modes and can be 
considered as a fixed classical background into which they exit [Hj. As discussed in 
^ in this limit the momentum quadrilateral constitutes a planar parallelogram (or a 
folded-kite). 
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Exactly the same reasoning applies where the exchanged particle is a scalar mode, 
and allows us to estimate the size of the contribution from such a process in the counter- 
collinear limit. Following Ref. [26] (see also Ref. [IS]), we write the background mode 
as ^^(ko). The effect of this mode is to re-scale the spatial coordinates so that the 
two-point function (CkiCk2) receives a correction proportional to its tilt, 

^CB(ko)(27r)35(ki + k2 - ko)Pdk2) ■ (3.1) 

A similar formula gives the transformation of (Ck3Ck4)- To obtain the contribution to 
the four-point function due to the scalar background mode (which we label with the 
superscript "SE") in the counter-coUinear limit, it follows that we need only to average 
over this pair of two-point functions, which yields 

((Ck.Ck2)cB(Ck3Ck,)cB)''^ = i2n)'S{J2K) ■ {ris - if ■ P^{h,) P^{h) P^{h) . (3.2) 

Thus, in the counter-coUinear limit, the correlations induced by scalar exchange are 
small, being suppressed by two powers of the scalar tilt. This is consistent with the 
estimate made in the Introduction (31]), which found the contribution from scalar 
exchange to be suppressed by 0(e^, ?7^, er^). Note that Eq. ( 13. 2p can be thought of 
as arising from the product of two back-to-back bispectra in the squeezed limit, and 
exactly satisfies the classical relationship between tnl and /nl described in Eq. (11. 6p . 
Nevertheless, the full trispectrum does not. Indeed, as we show below, although the 
contributions from the contact interactions are subdominant in this limit p8], those 
from gravitational wave exchange remain important and dominate the trispectrum. 

Now consider the effect of a long wavelength gravitational wave on (CkiCk2)- This 
will turn out to be much larger than the effect due to the scalar mode. In this case the 
effect of the background gravitational wave, 7^(ko), is to deform the geometry in such a 
way that a small spatial distance dx^ is rescaled to dx'^ {e'^)ijdx^dx^ ~ dx'^+'jijdx'^dxK 
In Fourier space we have 

(Ck.Ck2).. = (Ck.Ck2) + / ^7i^(ko) (^(Ck,Ck2)7.) , (3.3) 

where = k^ — ko/2. The effect of the background tensor mode is to deform the 
momenta such that k"^ k"^ — '-fijk^y [llj, so that the equation above becomes 

(CkiCk2)7B = (CkiCk2) - J ^^iiiko)klk] (^^^(c^^c^^)^^^ 

= (Ck,Ck2) - / |^(27r)35(ki + k,- ko)7i^(ko)^^Pc(^i) • (3.4) 

In this case, the effect of deforming the geometry is not suppressed by the scalar 
tilt, and therefore dominates compared with the effect of a background scalar mode. 
Nevertheless, Eq. (13.41) is proportional to the amplitude of gravitational waves, which is 
of order ~ r~^/^ when compared with the scalar amplitude. For this reason, one finds 
that the trispectrum in this limit is proportional to the ratio of the tensor to scalar 
power spectra. 
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Indeed, taking the average of (CkiCk2)7i3 and (Ck3Ck4)7s we find the relation 
((Ck,Ck.),,(Ck3Ck.),,)^^ = i27rmY.^a) 

a 

^ , Ui hi U urn q 

X E 4(ki2)^C(k34)^ ■ 4e ■ ^ ■ Pdh2)Pdki)Pdk3) • (3.5) 

5 1 3 

We can use Eq. f l2.32p to rewrite the polarization factor in this equation. In the limit 
fci2 0, ^1, 03 71 and Eq. fl2:32|) reads 

^etj{ku)k{k{eUku)kik^ = fc^A;^ cos 2x12,34 • (3.6) 

s 

Thus, Eq. (13.51) yields a semiclassical relation in the counter-collinear limit (fci2 -C fci ~ 
k2, ks ^ ki), 

(Ck.Ck.CkaCk,)^^ = (27r)35(J]kJ ■ ("Ar cos 2x12,34] ■ ^(^1)^(^3)^(^12) , (3-7) 

where we have used r = 16e and in this degenerate limit Xi2,34 becomes the angle 
between ki and k^. 

Let us check that Eq. (12.331) satisfies this relation and consistently reproduces 
Eq. (13. 7p in the counter-collinear limit. Indeed, for ki2 <^ /^i ~ k2,k3 ^ k^ the second 
line of Eq. (I2.33P dominates over the last two and we must only compute the sum of 
the integrals X1234 + X3412 by using Eq. ( 12.26^ in this limit. For ki2 = 0, ki = k2 and 
ks = k4. Using Eq. (I2.27p . 012 = 2ki, 034 = 2A;3, 612 = kf and 634 = k^. Thus, the right 
hand side of the first line of Eq. (I2.26P yields 2kik3, the second line yields kik^/A while 
the third line vanishes. Then, one finds X1234 +X3412 = ^kik^ and Eq. (12.330 becomes 

(Ck.Ck.Ck3Ck.)^^ = i^^mY.^'^) ^22Wk^ ■ ^ COS 2X12,34 . (3.8) 

a 1 3 12 

Comparing this equation with (13. 7p after using the relations P(;{k) = Hl/{Aek^) and 
r = 16e, one sees that the semiclassical relation is satisfied. 

Note that, strictly speaking, Eq. (13. 7p is not a consistency relation like Maldacena's 
consistency relation in the squeezed limit. In particular, the full trispectrum need not 
satisfy this relation. This relation will be satisfied by the full trispectrum only when 
contact interactions are negligible with respect to the graviton exchange in the counter- 
collinear limit. Accidentally, as we will show below, this is the case for slow-roll inflation, 
and the full trispectrum satisfies this relation. However, this is not the case, for instance, 
in models with non-standard kinetic terms. 

4. The total trispectrum 

The connected four-point correlation function of field perturbations at horizon crossing 
was computed in Ref. [28], assuming that inflation was driven by a collection of light 
scalar fields with canonically normalized kinetic terms. However, this calculation was 
restricted only to contact terms, neglecting correlations induced by particle exchange. 
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In this section we put together the trispectrum of Ref. for a single field with the 
exchange contribution computed in ^ This gives the full trispectrum to leading order 
in slow-roll. Then we discuss this result. 

When specialized to the case of a single field, the contact interaction (labelled with 
a superscript "CI") gives a contribution to the scalar trispectrum equal to [28] 

(Ck.Ck.CkaCk.)^^ = (2vr)^5(5^k„)^--|^ -^4(kl,k2,k3,k4) , (4.1) 

Q, LLa\ aJ pcfjjjs 

where we have defined a form-factor A^4 which encapsulates the detailed momentum 
dependence. This is given by the expression 

M,ik,, k2, k3, k4) = -2 ,}, I + 2k2 ■ Z34 + -^12^34 



"^12 "^34 "^t \ "^34 ^ 

where the quantities occurring here are defined by 

aab =ka-kb + kl , (4.3) 

Zab = CTabK - CTbakb , (4.4) 
TJ/ 1 I + kb 2kakb , . 

Wab = 1 + — ^ + , (4.5) 

jjr T , ka -\- kb -\- kc 2{kakb + feji^c + kgk^ Gk^kbkc , . 

Wabc = 1 + + ^2 + — ^ ■ (4-6) 

In order to write Eq. (14.11) from Ref. |28j, we have used the linear relation 
( = Lp^/\^ to express C iii terms of the field perturbations at Hubble crossing ^9^,. 
As discussed earlier, nonlinearities in the relation between ( and (f^, lead to corrections 
to the trispectrum that are suppressed by slow- roll with respect to (14. ip [28| HU [21]. 
Thus, at leading order in slow-roll, the total trispectrum is given by the sum of the two 
contributions, Eqs. (12.331) and (14. ip . 

(CkiCk2Ck3Ck4) = (CkiCk2Ck3Ck4)?^ + (CkiCk2Ck3Ck4)?^ • (4-7) 

We have already discussed the squeezed limit, in which one of the momenta becomes 
much smaller than the others. In this limit, the trispectrum must be described by 
Maldacena's consistency relation [HI |46l [26l [471 145]: the scalar mode associated with 
the small momentum forms a fixed classical background into which the remaining modes 
exit, in a similar fashion to the calculations discussed in ^ In Ref. [28] it was noted 
that the scalar contact term is one order lower in slow-roll than the term yielded by the 
consistency relation in this limit. These two apparently contradictory statements are 
nevertheless in agreement, because the leading-order slow- roll piece in the contact term 
vanishes in the squeezed limit. Although we have not pursued the calculation of the 
trispectrum to a sufficiently high order in the slow-roll to verify that Maldacena's relation 
holds in detail for the graviton exchange process, our final result for the trispectrum, 
Eq. (14. 7p . is also one order lower in slow-roll than predicted by the consistency relation, 
but vanishes in the squeezed limit. It is therefore consistent with Maldacena's result. 
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In the counter-coUinear limit ki2 <^ fci ~ A;2, ^3 ~ the form-factor Ai^ defined in 
Eq. (14.21) remains finite [28], so that the total trispectrum in Eq. (14.71) is dominated by 
the contribution from graviton exchange. One can show that in this limit it is possible to 
parameterize the total trispectrum (14.71) in terms of the local parameter t^l^'. Indeed, 
neglecting the contribution from contact interactions, the total trispectrum is given by 
Eq. (13.71) . Thus, comparing this equation with Eq. (11.91) gives 

4°r'=^^ cos 2X12,34. (4.8) 

Thus, in the counter-collinear limit, the trispectrum of single-field slow-roll inflation is of 
the local type in real space, with a magnitude expressed by Eq. (14.81) . This is consistent 
with our expectations: indeed, as was discussed in the Introduction (^, we expect 
classical processes operating outside the horizon to synthesize local contributions to the 
curvature perturbation. The counter-collinear limit of graviton exchange is tantamount 
to such a process: it corresponds to a gravitational wave which nucleates inside the 
horizon in the far past of inflation, and subsequently passes beyond the Hubble radius 
where it becomes approximately classical. In the diagrammatic interpretation, this 
extremely soft graviton then splits into four scalar quanta which propagate to their own 
hypersurface of horizon exit. These four quanta are correlated because they all originate 
from the same parent graviton, or gravitational wave. 

5. Discussion and conclusion 

We have calculated a new contribution to the inflationary trispectrum, arising from a 
process which can be visualized as graviton exchange. This contribution is of the same 
order of magnitude as the one previously calculated by considering only scalar and vector 
perturbations and gives rise to an effect of magnitude tnl ~ e- Although one might 
have expected t^l ~ 0{e'^,ri'^,eri) on the basis of Maldacena's consistency relation, we 
can conclude that the estimate tnl ~ e made in Ref. [28] is essentially correct, although 
modified in detail by numerical factors coming from the exchange process which was 
neglected in that work. Note that although we have framed our calculations in the 
context of single-field inflation, our results for the curvature perturbation also hold for 
models with multiple fields and a flat field space metric. 

We would like to emphasize that the relation tnl ~ /nl [see Eq. (II. 7p ]. which 
one might anticipate based on comparison of the local formulae (11.31) and (II. 6p . need 
not apply. Indeed, it does not apply in single-field slow-roll inflation models. As 
explained in the Introduction, one can only trust such classical relationships if the 
non-Gaussianity is generated by gravitational evolution subsequent to horizon crossing. 
In single-field inflation, non-Gaussianities are a quantum interference effect and are 
produced at Hubble crossing, after which the curvature perturbation ( remains constant. 

In Ref. [28j, it was supposed that the effect from quantum processes around horizon 
crossing would be maximized on a configuration close to the equilateral case, where all 
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ka take the same magnitude k. In the exactly equilateral limit \t^l\ can be computed 

analytically, giving a maximum effect which corresponds to 

23 r 
"^NL ^ 570^^'^ ~ 0.0231r ^ — (equilateral and g^^i^ = 0), (5.1) 

where to compute tnl from the trispectrum we have set ^^nl = 0. The tensor to scalar 
ratio r is currently subject to the observational constraint r < 0.22 at 95% CL from a 
combination of the WMAP five-year, baryon oscillation and supernova data [IS]. (This 
bound was approximated as tnl < r/50 in Ref. [28].) By searching through the available 
parameter space, we have subsequently located a marginally off-equilateral configuration 
which gives a non-linearity of greater magnitude, but opposite sign, 

T^l = 0.0298r ^ ^ (side ratios 1:1:1: 2.96 and ^nl = 0). (5.2) 

How large a correlation can be induced by the effect of graviton exchange? It is 
difficult to pursue analytic estimates for the graviton exchange contribution. Proceeding 
numerically, however, it is possible to find a marginally squeezed configuration which 
gives an effect corresponding to 

= -1.217r (side ratios 1: 10:10: 13.44 and ^NL = 0). (5.3) 

From Eqs. (15. 2p and (15.31) . it is clear that the total trispectrum will be dominated by 
the contribution from graviton exchange. Taking both effects together, we have located 
a configuration which gives a combined tnl equivalent to 

Tnl = -1.218r (side ratios 1 : 5.6 : 10 : 10 and ^nl = 0). (5.4) 

We would like to emphasize, however, that our numerical search has not been exhaustive, 
and it is possible that other configurations exist which give even larger ItnlI- These 
estimates hold for multiple-field inflation provided that the hierarchy among the sides 
of the momentum quadrilateral is not too large ^28] . For instance, the ratio of sides 
given in Eq. (15. 3p corresponds to an interval of ~ 2.6 e-folds between horizon exit of the 
modes corresponding to largest and smallest wavelengths. For such a small interval, we 
can expect that the approximation of almost simultaneous horizon exit is reasonably 
satisfied. On the other hand, in the limit where a large hierarchy exists among the ratio 
of sides in the quadrilateral, our analytical expression becomes untrustworthy. 

We can therefore conclude that the trispectrum at horizon crossing is generally 
dominated by the contribution from graviton exchange, and that at maximum the total 
trispectrum is roughly of order [tnlI ~ r. 
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Appendix A. Computing the time integral 



In this Appendix we compute the sum of the time integrals X1234 + X3412. Let us start 
by computing X1234, which reads 

Xi234 = r ^ r ^Im[(l + iA;i7y.)(l + iA;277*)(l-zM(l-^M)e'^''"''*^^'^+''^ 
J —00 V J —00 V 



X Im 



(1 + ik^'r]^){l + ikiri^){l - ikir]'){l - ik.A'q') 
X (1 + ik^2-n){l - iA;i27?')e*('''-''*)(*^^+'=^^e-^(''-'''^'=i2 



(A.l) 



It is convenient to divide this into the sum of 4 integrals. Using the rule Im(AS) 
-\{AB* - A*B)I2 we obtain 



X1234 — — ^ {I A — Ib — Ic + Id) , 



where 



n/2 



Ib 
Ic = 
Id = 



00 V J —00 V 



-00 V 
■00 

d7] 

-00 v"^ 



XY* 



-00 v'"^ 
" djf_ 

dr]' 
■00 V'^ 



Z* , 

z* 



(A.2) 

(A.3) 
(A.4) 
(A.5) 
(A.6) 



and 



(1 + ikir],){l + zk2V*){l - ikiv){l - ik2v)e'^'^~'^*^^''''^'''^ (A.7) 
(1 + zA';i2?7)e-^*% (A.8) 
Z{7]\7],) = {l+iks7j,){l+ik^7],){l - - - zA;i2V)e*^'''-''*)(*^^+*^^)+*''''i^ . 

(A.9) 

Each of these have to be solved separately. 

Let us first consider Ia- Expanding the T^^-dependent part near 77* — > 0, this can be 
written in the form 
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X r ^(1 -iair7 + a2^2 _i^^^3^)gi^ai 

J-oo T 

X f iPin' + P2r]'^ - if3sv")e''''^' , (A.IO) 

J —CO V 



where 



"l = ^1 + ^2 - ki2 , 


(A.ll) 


a2 = {h + k2)ki2 - kik2 , 


(A.12) 


"3 = ^1^2^12 , 


(A.13) 


/?1 = fcg + A;4 + ki2 , 


(A.14) 


/32 = -(/C3 + ^4)/Cl2 - ^^3^4 , 


(A.15) 


p3 = -^3^4^12 , 


(A.16) 


6 = ai + Pi . 


(A.17) 



The inner integral in the last line of Eq. (lA.lOP can be easily solved and its log divergent 
part cancels out. Plugging its solution into Eq. ( lA.lOp yields 



V a / 



( + ^ + - + i73 + 74^? + i75r?' ) e'^', (A.18) 



where the 7's are defined by 

P1P2 + P3 



7i = 


ttl — 






72 = 




73 = 


"3 - 




as 


74 = 






"3/^3 


75 = 





/3? 

^2 4 



(A.19) 



(A/32 + /93)-a2-^, (A.20) 
Pi 

(A.21) 
(A.22) 
(A.23) 



Carrying out the remaining rj integration, the final result is 



where i?i is the exponential integral or incomplete F function, 

r°° e"* 

Ei{z) = — dt. (A.25) 

Jz t 

This calculation can be repeated for the remaining integrals Ib,c,d- The final 
result (1A.24P can be generalized to Ib by letting ki —ki and k2 —k2', to Ic 
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by the replacements —k^, k^ —k^ and k^ —k^; and to lo by allowing all 

the ks to change sign. Eq. flA.2p can then be extended to X3412 by exchanging 12 34. 

Note that in the last line of Eq. ( 1A.24I) there are terms that diverge for r^^, — 0. 
However, we will now show that this line does not contribute to the final combination 
2^1234 + 2^3412 and can be thus dropped from the discussion. Indeed, the first and second 
terms in this line are the same for all the four integrals Ia,b,c,d- Thus they do not 
contribute to the sum in Eq. flA.2p . The third term of this line is odd under exchange 
of I A ^ Id and Ib ^ Ic- Thus, this term also does not contribute to the sum ( ]A.2I) . 
Furthermore, it is straightforward to show that the last term of this line does not 
contribute when considering the final sum X1234 + X3412. One might be concerned that 
Eq. flA.24p contains apparently divergent (and unphysical) terms in the limit 5 — > 0, 
which occurs when the momentum quadrilateral becomes counter-collinear. Although it 
is not possible to see cancellation at intermediate stages in the calculation, these terms 
also cancel in X1234 + X34i2 

We are therefore left with only the first line of Eq. (1A.24I) which can be computed 
straightforwardly, by considering the combined contribution from X1234 and X3412. This 
yields, using k^ = k^, 



X1234 + X' 



ki + k2 



3412 



a 



34 



-(034 + fci2)(a34 - 26, 



+ 



kik2 



634 
034 



ki2 + 



1-12 



ai2 



ksk^ — k 



34 



^34 



k^ 



12 



(^3 + k4 



12" 



034 



k 



12 



a^a^^kt 



612634 + 2kl2 n 



1 1 

kt ai2 
k 

+ 



+ (1,2 



+ (1,2 



3,4) 
3,4) 



12 



012034 ^t'3'12'3'34 / 



(A.27) 



where aab = ka + kb + kab and bat = {ka + kb)kab + kakb. The first line of Eq. ( ]A.27I) comes 
from summing all the contributions which correspond to the first term in the first line 
of Eq. ( ]A.24[) . while the last two lines come from the final three terms. 

* Note added after publication. We have chosen to nest the terms arising from (+, — ) and (— ,+) 
diagrams, which enter with a minus sign inside the square bracket of Eq. (|2.23p . The right hand side 
of Eq. (I2.23P could alternatively have been evaluated by factorization, which leads to an expression of 
the form 

((/:'k,^k,^k3V'kJ?^ = -J(2^)353(^k,)^e^(ki2)fclft^e,^„(k34)4fc4" 

a s 

2Ref r 5 r ^G<(^,,,y)G<J^.,,y)G<(^,,r;')G<(ry,,,?')^^fe<(^,'/) 

\^ —OC 'I J ~OQ 'I 

(A.26) 
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Once we explicitly replace the propagators in terms of mode functions, this expression yields no factors 
of S in the denominator. In this case, it is unnecessary to check that the limit (5 — > is regular and the 
calculation is simplified. We would like to thank Peter Adshead, Richard Easther and Eugene Lim for 
interesting correspondence on this issue. Note, however, that whichever method one chooses, the final 
answer is the same. 
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